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The author's lecture notes concerning the correlation functions and the thermodynamics of a simple polar 
fluid are summarized. The emphasis is on the dipolar hard sphere fluid and the mean spherical approximation 
and on the relation of these results to the Clausius-Mossotti and Onsager formulae for the dielectric constant. 
Previous excerpts from these lecture notes, Condens. Matter Phys., 2009, 12, 127; ibid., 2010, 13, 13002, 
have contained results that were not widely known. It is hoped that this third, and likely final, excerpt will prove 
equally helpful by gathering several results together and making these more widely available and recording a 
few new results. 
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1. Introduction 

This paper is dedicated to Yura Kalyuzhnyi on the occasion of his sixtieth birthday; it is the 
result of the beautiful work of Michael Wertheim on hard spheres and dipolar hard spheres that has 
inspired the author and many others, including Yura and his colleagues in Lviv. There is little in 
this paper, drawn from the author's lecture notes, that is not well-known to Michael but, perhaps, 
lesser mortals will benefit from this collection of results and those in two previous papers [TJ |2] 
taken from the author's lecture notes. 

Dipolar hard spheres are a simple representative molecular fluid and polar fluid. For a canonical 
ensemble, the well known definition for the h-particle correlation function for a simple atomistic 
fluid of N molecules in a volume V is easily generalized for a molecular fluid, 

V h Nl 



v n w i r 

9(1 ■ ■ ■ h) = Nh{N _ h)] ■ J exp[-/3C/(l • • • N)]dr h+1 ■ ■ ■ dr N dfl h+1 ■ ■ ■ dil N , (1) 



where 



N 



u(i---N)= u(Rii,toi,nj) (2) 
i<j=i 

is the energy of the system, (3 — 1/kT (T is the temperature and k is the Boltzmann constant), 
Qn is the configurational partition function given by 

Q N = J exp[-pU(l ■ ■ ■ N)]dr 1 ■ ■ ■ dr^dfii • • • dfl N , (3) 

and Rij — |r» — rv^| is the distance between the centers of a pair of molecules, i and j, whose centers 
are located at and Vj. The "volume" elements, dfl, are normalized so that J dfi = 1. Thus, 

d0= g inMg# 

47T 

The function u(Rij, Qi, Qj), that can be written as u(ij) for brevity, is the intermolecular potential 
between a pair of molecules. 
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The following notation is employed. A function, such as g(ij), that depends upon the orientation 
of each member of a pair of particles, i and j, is denoted by the presence of each of the indices of 
the two particles in the argument of the function. After integration over the orientations of the two 
particles, a function that depends only on the scalar separation, R, of the two molecules results. 
This spherically averaged function is denoted by the subscript 0. Thus, the radial distribution 
function (RDF) is given by 

9o(R) = y fl(12)dfiidn 2 = <(/(12)). (5) 
The Ornstein-Zernike (OZ) equation becomes 

h(12) = c(12) + p [ /i(13)c(23)dr 3 df7 3 , (6) 



where p — N/V, h(12) = g(12) — 1 and c(12) are the total and direct correlation functions. 
The common thermodynamic functions are given by 

E = ^NpJg(12)u(12)dr 2 dQ 1 n 2 = ^NpJ(g(12)u{12))dr2, (7) 

^ = 1 - <ff(12)«'(12)i2 12 )dr 2 , (8) 

and 



kT dp 



1 + p [ h (12)dr 2 . (9) 



dp 

In the above E is the energy in excess of the kinetic energy terms, p is the pressure, and /in(12) = 
<7o(12) — 1. The functions, h (12) and <7o(12), are the total and pair correlation functions, respec- 
tively. 

This article gives only an outline of the field of molecular fluids. The discussion will be restricted 
to molecular fluids with a hard core. For convenience, molecular fluids can be divided further into 
two broad types, (1) fluids in which the hard core is spherical (the asymmetry comes from the 
attractive potential) and (2) fluids in which even the hard core is nonspherical. The first class is 
conceptually simpler and is considered here. Dipolar hard spheres will be considered as an example 
of this class. Liquid crystals are an example of the second class and may, perhaps, be considered 
in a future installment. 

2. Dipolar hard spheres 

As an example of a molecular fluid, we consider the dipolar hard sphere fluid where the inter- 
molecular potential is given by 



CO, Rl2 < O", 

-jg-D{12), R 12 XJ, 



u(12) = { _> mio , (10) 



where p and a are the magnitude of the dipole moment and diameter of the dipolar hard spheres, 

£(12) = 3(ei • Ri 2 )(e 2 • Ri 2 ) - (ei • e 2 ), (11) 

where Ri 2 = Ri 2 /|Ri 2 | is a unit vector in the direction of Ri 2 , ei is a unit vector in the direction 
of dipole 1 , and 

A(12) = ei • e 2 . (12) 

The function A(12) does not appear in the intermolecular potential, except as part of D(12). 
However, -D(12) and A(12) contribute independently to the correlation functions. The dipoles are 
assumed to be nonpolarizable. 
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Barker [3] has proved the very useful theorem that is given in the following two equations, 

J (§i • a)dQi = 



and 



J • a)(e, ■ b)dQj = ~(a- b). 



Note that b could be e,, i ^ j. We shall call these results Barker's theorem. 

Using Barker's theorem, it is easy to show that 1, D(12), and A(12) are orthogonal, 



J L>(12)df7idfi 2 = 0, (13) 
J A(12)df7 1 d0 2 = 0, (14) 



and 



J L>(12)A(12)dOid^) 2 = 0. (15) 
The normalization of 1, D(12), and A(12) can also be obtained from Barker's theorem and is 

dQidQ 2 = 1, (16) 
D 2 (12)dSl 1 dQ 2 = ~, (17) 



and 



A 2 (12)dfiidft 2 = - . (18) 
3 



This means that 1, Z?(12), and A(12) are part of an orthogonal basis set. Indeed, they are a 
subset of the spherical harmonics. A basis set is a linearly independent set of functions with the 
property that any function can be expressed as a linear combination of the members of the basis 
set. One basis set for three-dimensional Euclidean space is the set of vectors in the directions of 
the x, y, and z axes. The space of all functions for which the spherical harmonics are the basis 
set has an infinite dimension. As will be seen, in some special cases the functions 1, £>(12), and 
A (12) form a complete basis set of finite (three) dimension but this is not usually the case. A 
basis set need not consist of orthogonal vectors or functions. However, it is convenient if they are 
orthogonal. A nonorthogonal basis set can be transformed into an orthogonal basis set by what is 
called the Schmidt orthogonalization procedure. Hence, we can expand 

5 (12) = g {R 12 ) + h D (R 12 )D(l2) + h A {R 12 )A(l2) + ■■■ , (19) 

where go(R\ 2 ) is given by equation (5), 

M-Ria) = \ J J D(12).g(12)dr! 1 d0 2 , (20) 
M fl ia) = 3 / A(12). 9 (12)d0 1 dr> 2 . (21) 



and 



The coefficients go(Ri2), /id(12), and h&{\2) can be called the "projections" of g(\2) onto the basis 
vectors, 1, D{\2). and A(12). The notation h, rather than g, is used for hp, and /ia because they 
are zero when R 12 — > oo. Note that h (R) = — 1, for R < a but ho{R) and h&{R) equal 0, for 
R < a. 
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The common thermodynamic functions are given by 

l + yo(<7)-^ 2 /^dR, (22) 
^JfeT-^/^dR, (23) 
and 

kT^ = 1 + p J h (R)dH, (24) 

where 2/(12) is the background, or cavity, function, y(12) = cxp[/3u(12)]<?(12) and u(12) is the pair 
interaction. Note that 2/(12) is a continuous function even if tt(12) is discontinuous. The functions 
yo(R) and ho(R) are the spherically averaged projections of 2/(12) and h(12), respectively. 

As we shall see shortly, the dielectric constant is also given by an integral involving h&. This 
means that the dielectric constant and common thermodynamic functions can be obtained from 
go, ho, and h\ even if the other projections are not known. Of course, in general, to obtain these 
three projections, the other projections must be calculated. In any case, these three projections can 
be called the active projections for the dipolar hard sphere fluid since they determine the common 
thermodynamic functions and the dielectric constant of this fluid. 



pV 
NkT 

E = 



3. Simple treatments of the dielectric constant 

The simple treatments considered here are based on the concept of the local field, E\ oc , felt by 
a dipole. This is not equal to the applied field, E, because of the other dipoles. Let us carve out a 
sphere of volume a, centered at a dipole. Since the dipole-dipole interaction is long ranged, we may 
assume that the dielectric or polar fluid is a continuum outside this sphere. Choose the volume of 
this sphere to be equal to the volume per dipole, 

4tt , V 1 

Y a -N-~ P - < 25 > 

The average value of p, is related to E\ oc , 

(m) - (pcosO), (26) 
/%os0exp(/3ju.E loc cos0)sin0d0 

W = ^t^ T^~v o\ ■ qao ■ ( 2? ) 

J cxp(PfiEi oc cos 0) sin 0d0 

As E\ oc is relatively small, the exponentials in equation (27) may be linearized. Thus, 

L w (l + yfyu.Ei oc cos0)cos0sin0d0 
J (1 + ppE\ oc cos 0) sin OaO 

The integrals of the first term in the numerator and the second term in the denominator vanish. 
The result is 

(m) = \Pn 2 E\ oc . (29) 
Our task is to calculate E\ oc . We will consider two simple approaches first. 

Clausius-Mossotti result for e 

The field inside the dielectric fluid is different from the applied field due to the polarization of 
this fluid. 

D = eE = E + 4ttP, (30) 
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Figure 1. Dielectric constant for the dipolar hard sphere fluid for pa 3 — 0.8 as a function of y. 
The solid curve gives the result of Tani et al. |23| . equation (142), and the broken curves give 
the CM Onsager [5], and MSA |12| results. The points are simulation results |5j- 



where D and P are the electric displacement and polarization vectors. For an isotropic system, the 
vectors have the same direction. Thus, 

P= e ^E = p{ l i). (31) 

Lorentz [4j (5j argued that there are four contributions to E\ oc : (1) the applied field, E; (2) the 
volume charge contribution of P, which is zero because P is a constant and V • P = 0; (3) the 
surface charge contribution of P on the surface of the sphere of radius a; and (4) the field due to 
the dipole, which is independent of E and so does not contribute to (/a) . For a surface element of 
this sphere at a polar angle, 9, measured from the direction of P and E, the area of this element is 
dS = 2na 2 sin 9d9. The surface charge density in the direction of the normal to the surface of the 
sphere at the polar angle 9 due to the polarization is P n — P cos 9. Thus, the element of the field 
due to the surface polarization is dE' = P n dS/a 2 and E' is 

E' = I cos0(Pcos0)27rsm0d0 = —P (32) 



so that 



4-7T 

E loc = E+—P (33) 



and 

(n\ = , . , . 

3^^ V 3 



(ri = lf3p 2 (E+^-P), (34) 



yielding 

P - 3_ 

Recalling that 

c — 1 J- - 

E ' 



p = ; 3 X7, Q .., - (35) 



e = l + ^, (36) 



this gives 

' ' 4 %/?M 2 = y (37) 



e + 2 9 

or 



l + 2y 
1-2/ 



(38) 
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which is the Clausius-Mossotti result [5j |7]. This is not a very good result because the CM e 
diverges when y = 1, for which there is no experimental support. Sometimes this problem is called 
the polarization catastrophe. The CM result for e is plotted and compared with some simulation 
results [8 in figure 1. There is no singularity in the simulation results. 

The constant y is not to be confused with the background function, j/(12). 



Onsager result for e 

To obtain E\ ocl Onsager [9] solved the boundary value problem for a sphere of radius a and 
dielectric constant equal to unity within an infinite dielectric medium whose dielectric constant is 
e and with an applied field E. Denote the potential inside and outside the sphere by fa and fa, 
respectively. Thus, 

V 2 fa = V 2 fa = 0. (39) 

The potential and displacement arc continuous across the surface of the sphere so that fa(a) = 
fa(a) and d(f>i(a)/dR — edfa(a)/dR. The potential fa is finite inside the sphere (in particular at 
R = 0) and, far from the sphere, fa = —ERcosO. The solution of this boundary value problem is 

fa (R) = - 7^j ER cos d - (4°) 



Hence. 



From this 



Using, 

the dielectric constant is given by 



(u) = -fiu 2 -^—E. (42) 
\w 3 pp 2e + 1 v ; 



^ = P(M> = (43) 



tyJh (44) 



or 



2e 

(e-l)(2e + l 



9e y - ^ 

Explicitly, e = [1 + 9y + 3y/l + 2y + 9j/ 2 ]/4. This is Onsager's formula. It is plotted in figure 1. 
This result does not diverge and is much better than the Clausius-Mossotti result. Until Wertheim's 
result, this was the standard formula. Wertheim's results will now be considered. However, some 
preliminary formulae are needed. 



4. Fourier transform of h^(R) 

As has been mentioned, h(12) can be expanded in spherical harmonics, 

h{\2) = h Q (R) + h A (R)A{12) + h D {R)D(l2) + ■■■ . (46) 

A similar expansion can be made for c(12). We will want to substitute these expressions into the 
OZ equation, equation (6). To do this it is convenient to use the Fourier transform. The Fourier 
transforms of ho, Cn, /ia> arid ca are straightforward. However, the Fourier transforms of hr> (12) and 
C£>(12) are more complicated because -D(12) contains R and we must transform the combinations 
h D {R)D{\2) and c D {R)D{\2) as wholes. 

First, recall that the Fourier transform pair is 

/(*) = / f(R) exp(ik • r)dr (47) 
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and 

f(R) = 7A3 / /(*) ex P(- ik ' r ) dk > ( 48 ) 

where R = |r| and k = |k| . Choose the coordinate system so that r = i?(sin 9 sin <j>, sin 9 cos <fi, cos 9) 
and k = (0,0, k). For an easier notation, define T(12) = ho(R)D(12). After some algebra, 

f = h D (k)D k (12), (49) 

where 

D k (12) = 3(ei • k)(e 2 • k) - (ei • e 2 ) (50) 

and 

00 

h D (k) = -4tt J R 2 j 2 (kR)h D (R)dR, (51) 


where 

, „ 3 sin x 3 cos x sin x 

,■ , • ^ ^ 

The function f(k) is sometimes called a Hankel transform. Note that hn{k) is not ho(k), the 
Fourier transform of hr>(R), which is given by 

00 00 
h D {k) = JT J Rs[nkRh D(R)dR = 4TT J R 2 j (kR)h D (R)dR. (53) 


The functions jo(x) and j2(x) are spherical Bessel functions. Equation (49) is a perfectly good result 
for the Fourier transform of hD{R)D{12) but it is a nuisance to have two kinds of transforms. Thus, 
it is useful to define an auxiliary function, 

00 

F(R) = f{R)-3 J ^PcLR', (54) 

R 

because, as is seen by straightforward integration, the Fourier transform of F(R) is the Hankel 
transform of f(R). 

Thus, in summary, the Fourier transform of T(12) = Hd(R)D (12) is given by 

f (12) = H D (k)D k (12), (55) 
where Dk (12) is given by equation (50) and 

00 

H D (R) - h D (R) - 3 J ^^IdR'. (56) 

R 

The inverse of the last equation is 

R 

h D (R) = H D (R) - A J H D (R')R' 2 dR'. (57) 


This auxiliary function has another interesting property. If f(R) is long ranged, F(R) is short 
ranged. For example, if 

f(R) = { °' x (5g) 

then 

^) = {o"^' nta. (59) 

This property can be exploited to evaluate integrals of long ranged functions which would be 
difficult if direct integration were attempted. 
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5. Fourier transform of the OZ equation for dipolar hard spheres 

It has been seen that 

h(12) = h (R) + h A (R)A(12) + h D {R)D{\2) + ■■■ (60) 



and 

Hence. 

and 



c(12) = cq(R) + c A (-R)A(12) + c D (R)D(12) + ■■■ . (61) 

h(12) = h {k) + h A (k)A(12) + H D {k)D k {\2) + ■■■ (62) 

5(12) = co(k) + 2 A (fc)A(12) + C D {k)D k {12) + ■■■ . (63) 

Table 1 . Wertheim's "multiplication" table. 





1 


A(23) 


D K (23) 


1 


1 








A(13) 





|A(12) 


±D K (12) 


D K (13) 





\D K {12) 


|[D K (12)+2A(12)] 



To take the transform of the convolution in the OZ relation, we must evaluate integrals of the 
form 

A(13)D(23)dft 3 . 

To do this, Wertheim's "multiplication table", which is given in table [T] is required. This multipli- 
cation table is easily obtained using Barker's theorem. Using this table, the transform of the OZ 
equation may be obtained. Since 1, A(12), and -D(12) are orthogonal, we can equate coefficients. 
Thus, 

h = co + ph c H , (64) 

h A = c A + -p(h A c A + 2H D C D ) + ■■■ , (65) 

and 

H D = C D + ^ P {h A C D + H D c A + H D C D ) + ■■■ , (66) 

with similar equations for the transforms of the higher order terms. 
We know that 

h {R) = -1, R < a, (67) 
h (R) -> 0, R -> oo, (68) 
Cn(-R) ->•(), R^oo, (69) 

h A (R)=0, R<a, (70) 
h A (R)~^0, R^oo, (71) 
c A (R)^0, R^oo, (72) 

h D (R)=0, R<a, (73) 
h D {R)^0, R^oo, (74) 

c D (R)^^-, i?^co. (75) 



and 
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Equations (73)-(75) are fine but we are interested in Hq{R) and Cd(R) rather than in /i£>(i?) 
and cd(R)- It is easy to show that 

H D (R) = -3K, R<a, (76) 
H D (R)^>-0, R^oo, (77) 
C D (R)^0, R^co, (78) 



where 



R 



The parameter, K. is independent of R but depends on T. p, e, etc., and is not known until the 
problem is solved. 

We can establish an interesting result for Cd(0)- We know that 

R 

c D (R) = C D {R) - J3 J C D {R')R I2 AR'. (80) 


Using equations (75) and (78) it follows that 

00 

J C D {R')R' 2 &R'. (81) 



R 3 R s 



From this, we have 



1 pC D (0)=y. (82) 



3 



6. Some exact results for e 

Onsager's expression is a special case of the exact result |10| ITT] 

(e- l)(2e + l) _ 4tt/3p(M 2 ) 
9e ~~ ~~ 9 N~' 

where M is the total dipole moment of the dielectric. Write this as 

(e-l)(2e+l) 



(83) 



9e V9K- (84) 



The parameter g^ is called the Kirkwood gx factor. The g^ factor can be written as an integral, 

9k = ^ - 1 + N{et ■ e 2 ) (85) 

yielding 

g K = 1 + J h A (R)dR = 1 + ^ph A (0) (86) 



so that 



(e-l)(2e+l) = 
9e y 



1 + ^a(0) 



(87) 



The Onsager approximation consists in neglecting the contribution of h A (R). 

Some other interesting exact results for e can be obtained using the OZ equation given above. 
We can use the truncated versions of the expressions for ho, h A , and hp, namely, 

h = c + ph c , 
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h A = c A + \p{h A c A + 2H D C D ), (89) 

and 

H D = C D + ^p(h A C D +H D c A + H D C D ). (90) 

The missing terms do not contribute. Solving the truncated equation (89) for h A gives 

c A + lpH D C D 
1 - \pC A 

Solving equation (90) for Hjj gives, for k = 0, 



h A = " A ^ 3 r . (91) 



Hd(0)= , d f\ 9 ^ 02) 
(x-y)(x + 2y) 



where x = 1 — ^/Oca(0) and equations (82) and (91) have been used. 
Equation (91) can be rewritten as 

1 r ^P 2 HdCd 
^ph A = 

3 1 - ^pCA 

Thus, 



i+yh A = f »_ u . (93) 



i + ^>Mo) 



1 - ^P^(0)y. (94) 



Using equation (92) yields 



X + V =l + ^A(0)= (£ - 1 f £ + 1) . (95) 



(x-y)(a; + 2y) 3 9ye 

The solution of this equation can be verified to be 



x = y e -±^ = l-^pc A (0) (96) 



or 

e- 1 y 



(97) 



e + 2 l-i p c A (0)' 

Hence, the Clausius-Mossotti result is obtained by neglecting c A (R). 
Finally 

= - ( ,-,K, + fr) (98) 

or 

^d(0) = - (£ ^. (99) 

These three routes to e may not be consistent for a given approximation. However, they will be 
consistent if the OZ relation is satisfied. 

Note also, that we have obtained exact expressions for h A (0), IJd(O), 5a(0), and Cd(0)! 

7. The mean spherical approximation for the dipolar hard sphere fluid 

Because the MSA is a linearized approximation, 1, A, and D are a complete basis set for the 
MSA. Thus, equations (88)-(90) can be employed. The MSA is 

(100) 
(101) 
(102) 
(103) 



ho 


= -1, 


R< 


CO 


- o, 


R> a 


h A 


= o, 


R < a 


c A 


- o, 


R> a 
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and 

H D = -3K, R<a, (104) 
C D = 0, R >cr. (105) 

These equations were first obtained and solved by Wertheim |12| . 

The first thing to note is that ho and Co are decoupled from equations (89) and (90). Equations 
(100) and (101) are the Percus-Yevick (PY) approximation for a hard sphere (HS) fluid. Thus. 

h (R) = h%l(R). (106) 

Algorithms for calculating h^g(R) = g^g(R) — 1, that are based on the formulae of Thiele [T3] and 
Wertheim |141 fTS] , have been given previously |16l ITT] . The other two equations may be solved by 
introducing the new functions, 



Hd{R) + \h A {R) 



(107) 



and 

h_(R) = ~[H D (R) - h A (R)]. (108) 
After a little algebra, the decoupled equations, 

h+=c+ + 2Kph+c+ (109) 

and 

h- = c--Kph-5- (110) 
follow. The MSA approximation consists of 

h+ = -1, R<a, (111) 

c+ = 0, R>cr (112) 

and 

h- = -1, R<a, (113) 

c_ = 0, R><7. (114) 

Hence, 

h+(R;p)=h^(R;2Kp) (115) 

and 

h.(R;p) = hlJ(R;-Kp), (116) 

where h^(R;p) are the PY hard sphere correlation functions. The equations for c + and c_ are 
similar. The contact values of h + (R] p) and h/R; p) are given by 

and 

h -^^ = - K "WT^' < 118 » 

where r\ = npa 3 /6. 

The algorithms of Smith et al. |161 fTT] can be used. These algorithms are robust and, with a 
small change, give sensible results, even for the negative densities required by equation (116). The 
required change is that the one line in the program where a cube root of a quantity involving rj 
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is taken, the instruction should be changed so that when rj is negative, the absolute value of r\ is 
used and the resulting cube root is multiplied by —1. 
Thus, 

h A (R) = 2K [h^(R; 2Kp) - h^(R: -Kp)] (119) 

and 

H D (R) = 2K hU(R;2Kp) + hH(R;-Kp) . (120) 



The contact values of h&(R) and Hd(R) follow from equations (117) and (118) together with 

i + v /2 



9ns <>; p) 



(i - 



(121) 



The function that we want is ho{R)-, not Hd{R). This can be calculated from equation (57). In 
particular, the contact value of ho{R) is 



h D (a) = H D {a)+3K. 



(122) 



The parameter K is not yet specified but we are in a position to do so now. Using equation (82) 



and 



yields 



1-P^(0) 



(1 + 2ry) 2 



(1-77)4 
(1 + 4Kr]) 2 (1 - 2Kri) 2 



3y, 



(123) 



(124) 



(1 - 2KT]) 4 (1 + Krj) 4 

which specifies K, which has been renormalized so that it is dimensionless. Note that < Krj < 1/2. 
When Krj = 0, y = and when Krj = 1/2, y = oo. 



o 



9j R > 




h A (R) 




1.0 



1.2 



1.0 1.2 

R/ff 



1.0 



1.2 



Figure 2. Correlation functions for the dipolar hard sphere fluid for pa 3 — 0.9. The points given 
by solid and open circles are the simulation results |18U19| for ftp? — (hard spheres) and j3p 2 = 
1, respectively. The solid and broken curves give the results of the MSA and equations (125) 
and (126), respectively, for /3p 2 — 1. 

The correction functions for the dipolar hard sphere fluid that follow from the MSA are plotted 
and compared with simulation results |181 HU] in figure 2 for a representative case. The MSA gives 
fairly accurate results for go(R). The simulation results for go(R) for dipolar hard spheres are very 
nearly equal to those for hard spheres but are slightly larger. Hence, one prediction of the MSA is 
that go(R) for dipolar hard spheres is independent of the magnitude of the dipole moment and is 
equal to the radial distribution function of a hard sphere fluid. This prediction is not exact but is 
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quite well satisfied by the simulation results. However, the MSA results for hr){R) and h&(R) are 
rather poor. Interestingly, the approximations, called LEXP, 

h D (R)=g^(R)h% SA (R) (125) 

and 

h A (R)=g^(R)hT A (R), (126) 

are much better. 



8. MSA thermodynamic functions 

Using the compressibility route, the thermodynamic functions are as follows: 

When the compressibility equation is used, this is a very poor result since the MSA incorrectly 
predicts, that there is no contribution from the dipolar part of the intermolecular potential. 
Using the pressure route, 

oo 

PV - 1 + 2 ^ 2 -^/^V (128) 



NkT {l-v) 2 3 J R 

which becomes 



PV -^ + ^ _ 3i ^ (129) 



and, using the energy route, 



NkT (1 - r]) 2 



E 3 

— r- = - - 3 K y- (13°) 

NkT 2 y y ' 




Figure 3. Free energy of the dipolar hard sphere fluid for pa 3 = 0.9 as a function of ftp 2 /a 3 . 
The points are simulation results [19) . The dashed curves marked MSA, 2, and 2+3 give the 
results of the MSA and perturbation theory when truncated after 2 and 3 terms, respectively. 
The solid curve gives the results of the Pade extrapolation of Rushbrooke et al. [24) . 



One interesting characteristic of the MSA is that the energy can be integrated analytically to 
give the (energy equation) free energy and this free energy can be differentiated to yield the (energy 
equation) pressure. After a little algebra, the results are 



A- A 



HS 



NkT 



= —Kt) 



(1 + Kr]) 2 (2 - Kr]) 2 



{l-2Kr j ) i (l + Kr)) 4 



(131) 
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and 



p-pm 
pkT 



= Kr] 2 



(1 + Kr]) 2 (2 - Krj) 2 



(l-2Kr]) 4 (l + Kr)) 4 



-3Ky, 



(132) 



where Ahs and Pus are the hard sphere free energy and pressure, respectively. Note that K has 
been renormalized so that it is dimcnsionless. The free energy that results from the MSA is plotted 
and compared with simulation results [20_ in figure 3. As is usually the case, the thermodynamics 
obtained from the energy equation are more accurate than those obtained from the compressibility 
or pressure equations. 



9. MSA dielectric constant 

The dielectric constant can be calculated by the three routes given above. All three routes yield 
the same expression for e. For example, starting with equation (97), we obtain 

e + 2 = (l + 4^ +2 (1-^) 2 , (133) 

e-l {\-2KriY {l + Krj) 4 y ' 

Solving for e gives 

(l+4i^) 2 (l + if77) 4 n _ 

6 = (i - 2KnY • (134) 

Some MSA results for e are plotted in figure 1. The agreement of the MSA result with the simulation 
results is better than for the CM and Onsager theories but the MSA results are still too small. 
Expanding the MSA expression for e gives 

£-,-£,» + ..., (135) 

which is correct to order y 3 . By contrast, the CM result is 

e - 1 



= V (136) 



and the Onsager result is, on expansion, 



1 



= V - W + • • • . (137) 

c T i 

There is no term of order y 3 in the CM theory. The Onsager coefficient of the y 3 term is too 
negative. 

Of course, approximations that are better than the MSA approximation can be used. For 
example, Fries and Patey [5T] used the hypernetted chain (HNC) approximation. The HNC results 
are better than the MSA results but the calculations are lengthy and, in contrast to the theories 
considered here, do not yield analytic results. It is to be noted that other combinations of /?, n, 
and p, besides y, appear when the HNC approximation is employed. This is true of other more 
general theories, for example the perturbation theory that is considered below. 



10. Perturbation theory for dipolar hard spheres 

Perturbation theory has been found to be very successful for simple fluids. It is natural to 
wonder if perturbation theory might also be useful for a polar fluid. The answer is yes but some 
qualifications are necessary. 

By expanding the free energy in powers of /3, the following result is obtained 

A = A HS + p 2 p 4 A 2 + p 3 p e A 3 + ■ ■ ■ , (138) 
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where Ajjs is the hard sphere free energy. The quantities A 2 and A 3 are given by 

P / 5r ffHSl^JdR = --p / Dfi dR (139) 



iVfcT A' J R 6 ynoy ' 6 ' J i? 6 

and 

^ 3 _ 1 2 /• (£>(12)D(13)£>(23)) , 100 ,, , _ ■ , ; 



v/ t ^ / <n p p 3 ffH S (123)dr 2 dr 3 - ^Iau , (140) 
TV/cT 6 7 (^12^13-^23)^ 54 

where 5hs(^) an d <?hs(123) are the pair and triplet distribution functions of the hard sphere fluid, 
and 

f 1 + 3 COS 0i COS 2 COS 3 , 10Q %, , / 1/f1 x 

^dd= / , p p p x 3 ffHs(123)dr2dr 3 , (141) 

where the 0, are the three interior angles of the triangle formed by the three sides, Rij. The angle 
0i is the angle opposite the side i?23j etc. Barker's theorem has been used to perform/simplify 
the angular integrations. The term of order j3 vanishes on angular integration, as do some of 
the terms of order f3 2 and /3 3 that formally contribute. Barker et al. [25] and Tani et al. [2"3~| 
have calculated I add by simulation and direct integration via the superposition approximation, 
<7hs(123) = 5hs(12).9hs(13)3hs(23). A numerical fit of their results is given by 



= 57r 2 j6 l + 1.12754p*+ 0.56192p* 

3 ° 1 - 0.05495/3* +0.13332p* 2 ' ' ~ 



where p* = pa 3 . 

As is seen in figure 3, this truncated series gives poor results. However the Pade, 

A = A ™ + i ,v l4h' (143) 

that was proposed by Rushbrooke et al. [24J, gives excellent agreement with the simulation results. 
A Pade tends to work best for alternating series. For example, the series 1 — 1+ 1 — 1+- • ■ is summed 
correctly to 1/2 by a Pade. Patey and Vallcau [20_ refer to the Pade results as "absurdly successful". 
This is meant as a positive comment and is a fair observation. Unfortunately, a Pade does not work 
well for the correlation functions. Some thoughts about the development of approximations that 
are consistent with equations (143) have been considered by Barker and Henderson [5S]. However, 
nothing much has come of these efforts. 

The dielectric constant can be calculated from 



(e-l)(2e + l) 

9e J 



1 QlddA 
167T 2 



(144) 



where 



f 3 cos 2 3 -l , 17tt 2 6 1 - 0.93952p* + 0.36714/9* 2 ,„ A . 

IddA = J jRMF 9m{U3)dr2dr3 = ^ l- 0.92398p*+0.23323p* 2 " (U5) 

This gives poor results, even with a Pade. However, the direct expansion, due to Tani et al. |23j, 



e = l+3 y + 3y 2 + 3 y 3 0^-l), (146) 

gives very good results, as is in figure 1. 

Perturbation theory can be recast by subtracting the MSA contributions from the perturbation 
terms and writing the perturbation theory as a series of corrections to the MSA. This was done by 
Henderson et al. [53]. The results are similar to those of the perturbation theory considered here. 
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1 1 . A few remarks 

The mean spherical approximation and perturbation theory are pleasing extensions of the 
classic theories of Clausius and Mossotti and Onsager for polar fluids. Not only do they make 
more accurate predictions for the dielectric constant of a polar fluid but they predict the other 
thermodynamic properties of polar fluids. Many of these ideas are applicable to polar fluids with 
a dispersion interaction. The simplest system of the kind is the dipolar Yukawa fluid. Szalai et 
al. [S7J and Mate et al. [25] have considered this model to be polar fluid. 

The molecules considered here are unpolarizable dipoles. Onsager considered polarizable dipoles. 
Valisko et al. have generalized some of the expressions presented here for polarizable dipoles and 
have made simulations for a polarizable dipolar hard sphere fluid. 

The author had hoped to include his lecture notes on liquid crystals as an example of a molecular 
fluid with a nonsphcrical hard core. However, despite some searching, these have not been found. 
If they do come to light, they can form a fourth part of this series and the "concerto" with three 
movements can become a "symphony" with four movements. 
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fle^Ki npocTi pe3y/ibTaTi/i p,j\n B/iacTMBOCTeR no/iapHMx rmi/miB 

fl,. TeHflepcoH 

Biflfliji xiMii' Ta 6ioxiMiT, yHiBepcuTeT BpireMa 3Hra, ripoBO, LUTaT \Ora 84602 

niflcyMOByeTbca /ieKL4iwHHW MaTepiaji aBTopa, npwcBaHeHwCi Kope/i3 L4i wh hm cpyHKL4inM i TepMOflMHaMiLii npo- 
CToro nojiapHoro njii/iHy. Oco6/inBa yBara npufli/iaeTbcn fluno/ibHOMy njii/iHy TBepflwx ccpep i cepeflHbo- 
ccpepuHHOMy Ha6.ni/i>KeHHK>, a t3ko>k 3B'a3Ky 14HX pe3ynbTaTiB i3 cpopMy.naMM K/iay3iyca-MocoTTi Ta 0H3arepa 
p/in flie/ieKTpwHHoT CTajioT. rionepeflHi Bi/iK/iaflKi/i i3 L4MX neKuji/i, Condens. Matter Phys., 2009, 12, 127; ibid., 
2010, 13, 13002, MicTMJin pe3ynbTa~™, aKi He 6y.n1/1 3ara/ibHOBiflOMi/iMH. € Haflis mo i\n TpeTH i, WMOBipHO, 
ocTaHHa Bi/iKJiaflKa 6yfle TaKoio >k Kopi/icHoto, o6'eflHyioHM KmbKa pe3ynbTa~riB i po6jiaHi/i fx flocTynHi/iMH fl/ifl 
Lui/ipujoT ayflWTopiT, a Ta ko>k npeflCTaB/intoHi/i KmbKa hobmx pe3y/ibTaTiB. 

K/iK>HOBi c/iosa: Kopenmj,itfHi cpyHKLtil] nonapm nntAHW, TepMop,tAH3M'mm dpyHKuJi, ftiejiexTpiAHHa CTana 
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